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ABSTRACT 


The  ground  plays  a  fundamental  role  in  the  study  of  the  ground  - 
vehicle-driver  complex  which  comprises  the  subject  of  .Land  Locomo¬ 
tion.  Through  the  pioneering  efforts  of  people  such  as  Col.  M.  G. 
Bekker,  formerly  of  the  Land  Locomotion  Laboratory  of  the  U.  S.  Army 
Ordnance  Corp,  Detroit  Arsenal,  great  progress  has  been  made  concern 
ing  the  role  of  ground  soil  properties  and  mechanics  in  determining  tin- 
tractive  ability  of  a  vehicle  traveling  on  the  surface.  However,  until 
quite  recently,  very  little  has  been  studied  concerning  the  role  of  the 
ground  surface  contour  in  determining  the  speeds  of  vehicles  especial¬ 
ly  on  off- road  conditions. 

The  authors  [3]  have  found  that  for  a  simplified  vehicle  traveling 
on  a  random  ground  contour  there  is  a  quantitative  relation  between  the 
ground  described  by  its  spectral  density  function  and  the  optimal  wheel¬ 
base  relative  to  the  criterion  of  maximizing  driver's  comfort.  This  was 
only  an  initial  study  but  points  the  way  to  the  fact  that  the  geometry  of 
the  ground  surface  must  be  clearly  studied  to  try  to  characterize  it  in 
such  a  way  that  quantitative  results  concerning  significant  vehicle  par¬ 
ameters  may  be  derived. 

.  -  1 

In  this  report  statistical  models  of  the  ground  surface  contour  are 
considered  along  with  possible  forms  for  the  two  dimensional  power 
spectral  density.  The  advantages  and  disadvantages  of  the  models  are 
pointed  out.  Perturbations  in  the  power  spectral  density  of  the  suriuce 


are  studied  to  determine  the  magnitude  of  their  effect  on  an  optimum 
criterion  chosen  for  vehicle  parameter  studies.  Finally,  the  power 
spectral  densities  from  an  actual  ground  surface  survey  are  presented 
and  discussed. 
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PREFACE 


This  is  the  fifth  in  a  series  of  papers  studying  the  motion  of  vehicles 
under  random  excitation.  Prior  papers  are: 

a.  Land  Locomotion  Report  No.  48,  Behavior  of  a  Linear  One  Degree  of 
Freedom  Vehicle  Moving  With  Constant  Velocity  on  a  Stationary 
Gaussian  Random  Track. 

b.  Land  Locomotion  Report  No.  56,  On  the  Behavior  of  a  Linear  Two 
Degree  of  Freedom  Vehicle  Moving  With  Constant  Velocity  on  a 
Track  Whose  Contour  i3  a  Stationary  Random  Process. 

c.  Land  Locomotion  Report  No.  65,  On  the  Statistical  Analysis  of  the 
Motion  of  Some  Simple  Vehicles  Moving  on  a  Random  Track,  which  also 
includes  an  appendix  which  defines  many  terms  and  derives  many 
basic  concepts  used  in  the  series. 

d.  Land  Locomotion  Report  No.  66,  On  the  Statistical  Analysis  of  the 
Motion  of  Some  Simple  Two-Dimensional  Linear  Vehicles  Moving  on 

a  Random  Track. 
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I.  Introduction 


The  remarkable  speed  range  of  modern  passenger  and  commer¬ 
cial  motor  vehicles  on  a  modern  highway  is  due  partly  to  the  power 
available  in  contemporary  power  packages,  partly  to  the  suspension 
system,  partly  to  the  straightness  of  the  highway,  but  it  is  mainly 
due  to  the  smoothness  of  the  surface  upon  which  they  operate.  One 
immediately  realizes  the  importance  of  the  smoothness  of  the  surface 
when  these  vehicles  are  required  to  operate  off  the  road.  Certain 
wheeled  military  vehicles  have  nearly  equivalent  speed  ranges  when 
operating  on  the  same  types  of  smooth  pavements.  It  must  be  remem¬ 
bered,  however,  that  of  necessity,  a  military  vehicle  may  be  required 
to  operate  off  paved  roads.  Under  off-road  conditions,  these  military 
vehicles  like  modern  passenger  and  commercial  vehicles  have  a  very 
limited  speed  range.  This  fact  is,  of  course,  not  strange,  since  the 
same  considerations  are  used  in  the  design  of  the  three  types  of  vehicles. 

The  basic  design  criteria  for  all  wheeled  vehicles  in  most  categor¬ 
ies  is  much  the  same--good  speed  and  control  on  the  highway  together 
with,  in  some  cases,  load  carrying  capacity.  A  vehicle  designed  using 
such  criteria  does  not  operate  and  cannot  be  expected  to  operate  at  even 
moderate  speeds  when  the  assumption  that  it  will  operate  on  a  modern 
smooth  highway  is  changed  to  the  assumption  of  operation  under  off¬ 
road  conditions. 

To  secure  reasonable  off-road  speeds  of  wheeled  or  other  types 
of  vehicles,  it  is  obvious  that  fundamental  studies  must  be  made  of  the 
factors  which  limit  speed  under  these  conditions. 
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The  Land  Locomotion  Laboratory,  Detroit  Arsenal,  has  ad¬ 
vocated  that  such  studies  be  made.  The  material  presented  in  this 
report  is  results  of  such  studies  conducted  by  the  Midwest  Applied 
Science  Corporation  of  Lafayette,  Indiana. 

This  report  in  particular  is  concerned  with  the  statistical 
description  of  the  ground  surface,  and  possible  methods  of  char¬ 
acterization. 

We  hasten  to  add  that  this  is  a  preliminary  study  of  the  prob¬ 
lem.  Therefore,  the  results  must  be  considered  as  tentative  and 
apt  to  change  as  more  terrain  is  surveyed  and  the  data  is  analyzed. 


II.  The  Ground  -  Introductory  Remarks 


It  is  a  remarkable  fact  that  although  the  very  earliest  form  of 
transportation  was  land  locomotion,  the  way  in  which  the  ground  char¬ 
acteristics  should  enter  into  the  design  of  vehicles  has  been  the  least 
understood  of  the  three  major  environments  of  travel,  that  is,  of  land, 
sea,  and  air  travel. 

We  shall  be  primarily  concerned  here  with  the  ground  part  of 
the  ground-vehicle  driver  complex.  Any  study  of  the  Theory  of  Land 
Locomotion  must  start  with  the  ground  for  it  is  the  primary  and  con¬ 
tinual  input  into  any  vehicle  that  moves  by  way  of  its  contact  with  the 
ground,  whether  the  contact  is  through  tires,  tracks,  skis,  or  any 
other  mode. 

The  characteristics  of  the  ground  can  be  broken  up  into  two 
major  areas,  the  physical  properties  of  the  ground  as  a  material  and 
the  geometrical  properties  of  the  ground  as  a  two  dimensional  surface. 

The  material  properties  of  the  ground  are  those  associated  with 
its  elasticity,  plasticity  and  tr affic ability .  Bekker  £1J  *  and  his  as¬ 

sociates  have  made  many  interesting  and  important  studies  of  these 
ground  properties  and  their  relation  to  the  Theory  of  Land  Locomo¬ 
tion. 

It  is  the  ground  as  a  geometrical  surface  that  has  been  the 
least  studied  relative  to  the  Theory  of  Land  Locomotion.  Furthermore, 


the  surface  apparently  has  not  been  analytically  characterized  in 
any  realistic  fashion.  This  can  be  thought  to  be  mainly  due  to  its 
irregular  nature  and  the  lack,  until  relatively  recently,  of  methods 
of  analyzing  and  characterizing  surfaces  with  such  peculiarities  or 
--if  you  will- -randomness  in  elevation.  The  methods  we  refer  to 
are  those  found  within  the  scope  of  the  Theory  of  Probability  and  Ran¬ 
dom  Functions . 

Only  recently  have  statistical  analyses  of  road  surfaces  been 
considered,  in  particular  the  work  of  Hoboult  et  al  [  2]  in  determin¬ 
ing  the  spectral  properties  of  airport  runways.  Even  here  ground 
vehicles  were  not  the  primary  consideration! 

In  the  studies  over  the  past 'two  years  by  Bogdanoff  and  Kozin 
of  the  Midwest  Applied  Science  Corporation  [  3]  [4]  under  con¬ 
tract  with  the  Land  Locomotion  Laboratory  of  the  Detroit  Arsenal 
OTAC,  U.S.A.,  ground  elevation  surveys  have  been  performed  along 
lines  on  typical  plowed  ground  in  central  Indiana.  Spectral  analyses 
were  then  performed  from  these  surveys.  Within  their  studies  the 
statistical  analysis  was  performed  of  simple  vehicles  traveling  on  a 
track  whose  contour  is  one  member  function  of  a  random  process  hav¬ 
ing  spectral  properties  approximating  those  found  from  the  surveys. 

Full  details  of  these  studies  are  given  in  the  next  paper.  At  the  time 
of  the  writing  of  this  paper  the  results  of  these  surseys  coin  <  rnir.g  » In¬ 
field  as  random  surface  have  not  been  fully  appraised.  Therefore, 
even  though  the  statistical  properties  are  known  along  partn  alar  st  r.ogir 
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lines  on  this  field,  its  general  two-dimensional  statistical  character 
remains  unknown.  It  is  probably  safe  to  say  that  the  question  of  the 
statistical  character  of  arbitrary  ground  surfaces  is  quite  new  and 
yet  to  be  answered. 

In  this  paper  we  shall  consider  two  aspects  of  the  theory  of 
random  ground  surfaces,  the  theoretical  and  the  experimental.  Un¬ 
der  the  theoretical  considerations  we  shall  discuss  a  possible  anal¬ 
ytical  model  of  the  random  ground  surfaces  along  with  the  associated 
two  dimension-spectral  density  functions.  The  advantages  and  short¬ 
comings  of  using  these  models  will  also  be  discussed.  It  is  certain 
that  approximations  to  ground  surfaces  through  spectral  densities 
will  have  to  be  employed  in  studying  the  motion  of  vehicles  traveling 
over  random  surfaces.  Hence,  the  magnitude  and  nature  of  the  errors 
introduced  into  the  analysis  by  using  such  approximations  will  also  be 
studied. 

Under  the  experimental  considerations  we  shall  present  and  dis¬ 
cuss  the  results  of  a  ground  surface  survey  that  was  designed  specific¬ 
ally  to  try  to  characterize  the  spectral  properties  along  straight  linen 
in  any  direction  for  that  particular  field.  This  survey  will  be  compared 
with  the  theoretical  results  discussed  in  this  paper  in  order  to  determ¬ 
ine  the  feasibility  of  a  spectral  characterization  of  the  type  we  propose. 
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Ill.  The  Ground  -  Random  Models 


The  model  of  a  random  surface  that  has  gained  great  favor 
among  oceanographers  and  meteorologists  studying  the  surface  of  the 
sea[5]  [6]  is  simply  the  two  dimensional  analogue  of  the  Rice 
random  functions  [7}  that  are  used  as  models  for  the  random  tracks 
in  the  paper  to  follow.  To  this  end,  consider  the  expression 

D  z  ( x,y)  m  £  ck  exp  i  (uk  x  +  vk  y+  <pk)  , 

km-n 

where  ck  *  uk  *  vk  are  rea*  numbers  satisfying 


vk  ,  uk  ,  Vjf  hove  the  dimension  L 

The  <p{s  are  independent  random  phases,  each  uniformly  distributed 

on  the  interval  (  0,2w)  satisfying  — ^  .  Hence,  we  see  that 

x  (x,y)  is  real  and,  indeed,  is  no  more  than  a  linear  combination 
of  cosine  terms  with  random  phases.  Conceptually  this  fits  the  pic¬ 
ture  that  one  has  of  a  ground  surface  full  of  miniature  mounds  and  val¬ 
leys  . 

One  immediate  objection  to  the  form  1),  is  that  as  n  — *•<»  , 

Z  ( x,y)  must  become  a  Gaussian  Process  [6]  However,  it  is  strong¬ 
ly  suspected  that  the  ground  surface  in  general  is  not  Gaussian. 

It  is  known  that  except  for  a  very  small  error  due  to  non-linearities 
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in  the  hydrodynamic  equations  [5]  the  sea  Surface  is  Gaussian 
in  character  and  hence  1)  modified  to  include  time  variation  fits 
the  phenomenon  quite  well, 

but  this  objection  can  be  dispensed  with  since  our  studies 
f  eh  Y  dynamics  have  beer  directed  toward  analysing  the  so- 
>  ailed  :-ei  ond-orrier  properties.  These  are  properties  that  are 

r,m ”i  d  by  the  cov  ariance  function.  It  can  easily  be  shown.  [81 
that  when  studying  second  orde--  properties  of  an  arbitrary  second 
order  random  process,  that  is  a  process  whose  mean  square  exists, 
the  process  may  be  replaced  by  a  Gaussian  process  with  the  same 
covariance  function.  Thus  no  contradiction  is  created  by  using  1). 
We  realize  fully  that  l)  may  not  suffice  if  we  must  study  Subtle 
properties  of  the  vehicle  motion  that  are  based  upon  considerations  of 
the  probability  densities  themselves  and  not  second  moments. 

Upon  taking  the  first  two  moments  of  z  ( x,yJas  given  by  1  ) 
we  find  r  1 

E  [ z(x,y)\  =  0 

3)  E  [ z(x,y)  z  (x  +  x1,  y+y')  1  c£expi(uk  x'+vk  y')  =f(xi  yi) 

•*  k=-n 

where  EC  3  denotes  the  expectation  or  ensemble  average  of  the 
quantity  within  the  brackets. 

In  words  3)  states  that  the  mean  is  zero  and  the  covariance 
function  exists,  being  dependent  only  upon  the  differences  of  the  ab- 
c  issas  and  ordinates.  Therefore,  by  definition,  z(x,y)  is  second 
order  and  at  leasi  ,\eakly  stationary  random  process. 
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Because  of  the  weak  starionarily  expressed  by  3),  one  may 


define  its  Fourier  inverse  as 

4  F(  x,y  )  -  jdu ^dv  %  (u,v)  exp  i  (ux  +  vy) 


(Note:  Unless  explicitly  stated  otherwise,  all  integrals  will  be  over 
the  range  (-oo,ao)) 

Hence, 

5)  %  (uv)  =  Jdx  Jdy  r ( x,y)  exp  i(xu+yv) 


The  function  t(u,v)  is  the  two-dimensional  power  spectral 
density  function.  That  is  f(u,v)Au  dv  yields  approximately  the  aver¬ 
age  amount  of  power  of  the  random  process  that  is  found  in  the  fre¬ 
quency  rectangle  ( u,u  +  A  u)  X  (v,  V  +  &v  ) . 

We  may  write  explicitly, 


6) 


where  the  index  k'  runs 


over  those  frequency  pairs  ( )  that  are  in  the  rectangle 

(u,u  +  6u)  X  (v,v+av). 


This  :s  mi!  i  i\  tin-  i  intension  to  two  dunen.s  ions  of  the  power 
spectral  density  tunition  that  one  encounters  when  studying  one-di- 
mensionaJ  records. 

A  nmjor  question  may  now  be  asked.  The  question  is,  ‘'How 
does  one  determine  a  useful  approximation,  economically,  in  both 
time  and  expense,  of  £(utv)  for  an  actual  ground  surface?" 

Furthermore,  what  can  be  said  about  t  (u.v)  restricted  to 
a  given  path  on  the  surface? 

The  reason  we  feel  that  these  are  major  questions  is  due  to 
the  tact  that  the  power  spectrum  of  the  surface  appears  to  be  the  most 
readily  available  way  of  characterizing  the  surface  statistically. 
Whether  or  not  one  may  be  able  to  determine  simpler  characteristics 
of  a  random  surface  that  can  be  used  significantly  for  vehicle  design 
purposes  seems  to  be  unknown  at  present.  Much  research  remains  to 
be  performed  in  this  area. 

The  remainder  of  this  paper  will  be  devoted  to  trying  to  answer 
t  tie  questions  posed  abo^e  . 
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i  V  .  Spectral  PrupfrlU’S  Along  Paths 

We  shall  first  consider  straight  line  paths  through  arbitrary 
points  ( xQ,y0)  having  directions  determined  by  the  angles  formed 
with  the  positive  x  -axis,  as  shown  in  Figure  1. 


Figure  1. 


7) 


Following  the  procedure  in  [6]  consider  the  transformation 

£  =  ( x—  x0) cos  0+( y  -  yQ  )sin  0 

T ; 

V  =-(x  -  xQ)  sin0+(y  -yQ)cos  0  , 


and  its  inverse 


X  =  x0  +  £  cos  0-  7}  sin0 
y  =  yo  +  £  sin  0+  7}  cos  0 


The  line  in  question  is  given  by  *f  identically  equal  to  zero. 
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Let  us  first  notice  that  for  a  weakly  stationary  protest  along 
a  line,  l  ,  we  have 

8)  rnr>)=  E  [z(T)  2  (T+Y')\ 

where  V  —  V-  k  U  and  /  U  /  —  ft  *£,  *  U  ~  0  • 

This  is  given  schematically  in  Figure  2. 


Figure  2. 


Hence  _  .  . 

9)  V  =  k  U 

Thus  we  see  that  the  covariance  along  any  line  depends  only  upon  the 
angle  S  ,  or  the  direction  of  the  line.  Therefore  all  classes  of  paral¬ 
lel  lines  have  the  same  covariance  function,  and  it  follows  that  they 
shall  have  the  same  spectral  density  function. 
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No'.v,  using  (4  )  and  (7)  withf*^,  y0)—(0,0), we  find 


r'(Cv>  = Jdu Jdv  •  &(u,  v)  exp  /£  (  ( cos&-  7  sin  &Mfsin&+  7?cos&)v j  j 
=■  Jdu Jdv {j£ ( u, v)expi  cos&+  v sin&)  +  7}(—usin&+vcos  <*»]} 

10)  -  =  JduJdv'  (u‘t  v')  expl  (£ u'+tjv')  ,  where 

u' 3=  u  cos  &+  v  sin  & 
v'=—u  sin  0+v  costf 

^  $  (u  (u'.v') ,  vfu'y)),  and  the 

Jacobian  of  the  transformation  is  unity. 

Placing  jm  0  for  the  line  in  question  yields 

a)  IJ  (£)  =  r  ({,0)  =  Jdu'Jdv' Ja  (u't  v')  exp  i  £ u' . 

But  the  covariance  and  spectral  density  functions  are  Fourier  trans¬ 
form  pairs,  so  that, 

12)  %Ju')  =  Jdv'j* '  (u'.v')  -  Jdv'  $(u,v)  ,  which  is 


the  spectral  density  function  along  the  line  in  the.  given  direction. 

At  this  point' we  wish  to  strongly  stress  the  fact  that  all  of 
the  results  so  far  obtained  in  this  section  are  based  only  upon  the 
weak  stationarity  of  the  randcm  surface  process  and  not  upon  explicit 
forms  or  processes. 
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From  11),  we  obtain  an  interesting  result  concern. ng  weakly 
stationary  random  surfaces. 

i3j  (O)  =  fdu'fdi/'  %!(u',  v‘) 

-  Jdu  Jdv  £(u  ,v) 

But  the  last  double  integral  is  independent  of  0  .  Therefore,  the  var¬ 
iance  must  be  the  same  along  any  direction!  Apparently  this  has  been 
found  to  be  the  case  for  sea  surfaces,  but  it  is  not  clear  that  one  shall 
find  this  to  be  true  for  the  ground  surface.  The  consideration  of  weakly 
stationary  random  processes  to  describe  the  ground  surface  could 
break  down  seriously  at  this  point. 

Let  us  now  study  paths  that  are  not  necessarily  straight  lines. 
What  can  be  said  about  their  power  spectral  density  functions  ?  We 
shall  approach  the  question  by  first  considering  the  covariance  func¬ 
tion  along  a  path. 

We  shall  parametrize  the  path  using  arc  length,  s,  along  the 
path  as  the  parameter.  The  elevation  for  the  path  (s) ,  y(sj)  on  the 
surface  1)  is 

t4)  z  (x(s),  y(5))  =  Z  ck  e*pi\uk  x(s)i-\  y(s>  +  <Pk  } 

k=~n 

The  mean  is  still  zero. 
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For  the  covariance,  one  finds, 


£  \z(x<s) ,  y(s))  2  (x(s+s'J  ,  y  (s+s'J)^ 

-  £  [  £  £  £  Cj  expi(uk  x(s)+vk  y(sj+uj  xfs+s'J+Vj  y(s+s')+<f>.  +<Q) 

k—n  J—n 

—  £  cfexpi^Uj  [jrfstsy- x(s)^+vj  fs+s‘)-y(s)  J  } 

Jm-n 


since  there  can  only  be  a  contribution  to  the  average  if  k»  —  j 
But  in  order'for  this  to  be  a  function  of  S'  alone  so  that  a  spectral 
density  function  exists  along  the  path,  it  follows  that 

x  ( s+s')-  x  (s)  —  g(s') 

/  ft )  * 

y  (s  +  s')~  y(s)  =  k(s') 


for  any  s ,  s'  along  the  path. 

The  obvious  example  of  16)  is  the  straight  line  path  that  we 
have  already  considered.  Are  there  any  other  cases?  The  answer 
is  easily  seen  to  be  no  by  simply  differentiating  the  relations  16)  with 
respect  to  s  and  setting  s  equal  to  zero,  yielding 

17)  dt  '  ~<fcS  ccns,an>- 

Therefore,  a  spectral  density  function  cannot  he  expected  to  dura^  - 
teri/.e  paths  other  than  straight  Ju.es  on  a  weakly  stationary  random 
surfac  e  . 
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V.  Spectral  Models  for  Real  Ground  Surfaces 

From  the  preceding  considerations,  we  know  that  spectral 
densities  exist  along  straight  lines  on  weakly  stationary  random 
surfaces.  We  may  assume  that  fields  free  from  streams,  shell  holes 
or  other  man  made  singularities  such  as  irrigation  channels,  etc. 
have  the  required  stationary  properties,  so  that  spectral  densities  of 
the  ground  elevations  along  straight  line  paths  can  be  obtained.  Our 
problem  in  this  section  is  inverse  to  that  considered  in  the  previous 
section,  namely,  how  can  we  determine )l(u,v)  from  Spectral  proper¬ 
ties  of  straight  line  paths  ? 

From  several  samples  of  surveys  taken  along  actual  rough 
ground,  it  appears  that  the  power  spectral  density  function  has  the 
general  form  for  straight  lines  as  given  by  the  exponential  function 

/8)  1$T  •**  (~  ~£f)  ■ 

u0  is  a  measure  of  the  spread  of  frequencies  about  zero  and 
<T *  is  a  measure  of  the  area  under  the  curve,  that  is,  the  variance. 
We  hasten  to  add  that  18)  is  by  no  means  definitive,  since  the  statis¬ 
tical  characterization  of  ground  surfaces  is  still  in  its  infancy.  Qual¬ 
itative  results  may  be  based  upon  the  form  18),  but  quantitative  re¬ 
sults  based  upon  18),  must  contain  inaccuracies.  The  magnitude 
and  effect  of  errors  based  upon  inaccuracies  in  the  estimate  of  the 
power  spectral  density  function  will  be  considered  in  the  following 
section. 
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Upon  initial  considerations,  there  are  two  possible  ways  of 
extending  the  function  18)  into  a  two  dimensional  spectral  density 
function.  Each  of  the  extensions  represent  distinct  inodes  of  attack 
on  the  problem  of  obtaining  the  )$(u,v) function  in  the  easiest  possible 
fashion.  The  functions  extensions  are 


19) 


*f,(u,v)-  —  YfP-zJiL-  0Xp-  ^au'+Zbuv+cv*} 

V"'  -  uj&fhr  •*'  ~  T^ffr 


The  first  expression  is  simply  an  extension  of  the  one  dimen¬ 
sional  exponential  function  to  a  two  dimensional  one.  The  second  ex¬ 
pression  represents  the  variation  in  18)  along  linear  paths  defined 
by  9  the  angle  between  the  line  and  the  positive  n  -axis. 

To  explicitly  determine  the  functions  19)  one  would  need 
a ,  b,  c  in  the  first  case.  In  the  second  case  the  variation  of  uQ  with 
9  would  be  needed.  The  factor  is  an  amplification  factor  that 
represents  the  variance.  It  must  be  constant  along  any  line  as  deriv¬ 
ed  above.  Notice  that  t  /V, ^-constant  in  the  first  expression  of  19) 
correspond  to  ellipses  in  the  (u,  v  ) -plane. 

To  help  substantiate  the  form  chosen  for  the  first  of  19),  we 
notice  that  one  may  define  moments  for  t(u,v)  as  £63 

20)  rnpq  m  Jdu  Jdv  $(u,v)  up vq  . 
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This  allows  us  to  define  moments  along  straight  lines  in  various 
directions  by  making  use  of  the  transformation  given  in  1.0). 

2i)  mn(&)  =  fdo1  Jfy(u')  u,n 

=  jdu' fdv‘ Z  fu'.v') u,n 
=  Jdu  JdvJZ  (u,  v)  (u  cos  &+v  sin &)n 

Notice,  in  particular,  that  the  second  moment  is 

22)  rng  (&)  =  mgQ  cos*&+  2  mn  cos  &  sin  &+  mQg  sin*&. 

This  function  of  9  is  continuous  in  its  domain  of  definition 
and  must  attain  its  extreme  values.  However,  determining  these 
values  is  equivalent  to  finding  the  stationary  values  of  the  quadratic 
form 

23)  m*o  x*+2muxy+m0gy 

subject  to  the  conditions  x  +  y  —  / 

The  solution  to  this  classical  eigen- value  problem  is  given  by 

mto-  A  tnn 

24) 

mn  moi~  * 

We  have,  then,  two  distinct  directions  for  which  mg(&)  is  sta¬ 
tionary.  These  directions  can  be  considered  to  correspond  to  the  di¬ 
rections  of  the  major  and  minor  axes  of  the  elliptical  power  spectral 
density  .distribution. 
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To  determine  the  elliptical  power  spectral  density  function  for 
a  givenfieJd  One  would  have  to  (Ificrnunr  a ,  b  ,c  .  Hhis  can  be  accom¬ 
plished  by  surveying  three  linear  paths  along  three  distan  t  directions. 
We'  shall  discuss  a  survey  devoted  to  this  task  later. 

The  second  of  the  forms  in  19)  can  be  useful  if  there  is  a  def¬ 
initely  preferred  direction  of  the  variations  in  contour  of  a  particular 
ground  surface.  This  corresponds  to  the  so-called  long  crested  prop¬ 
erty  in  a  wave  pattern.  Such  a  state  would  presumably  exist  shortly 
after  a  field  had  been  plowed,  or  if  there  happens  to  be  a  prevailing 
weather  and  wind  direction  creating  a  distinct  grain  to  the  surface. 
Figure  3  represents  schematically  this  effect. 
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Therefore  for  a  wave  length  of  l  in  the  preferred  direction, 
we  find  that  the  length  becomes  — -yp  along  the  direction  defined 
by  0,  0  <  -j-  • 

The  original  frequency  u  becomes  u  eos0  along  the  new  di¬ 
rection.  If  the  power  spectral  density  function  along  the  preferred 
direction  is  t  (u)  then 


Hence,  using  25)  the  second  of  19)  must  become 


26) 


uQ  cos&\/2w  eXf>  2  u*cos*& 


Therefore  becomes  uQcos  0.  We  now  have  the  spectral  prop¬ 

erties  along  any  linear  path  on  the  surface  merely  from  the  results 
found  along  the  preferred  direction.  For  fields  that  can  be  recognized 
to  have  the  properties  mentioned  above  only  one  survey  shall  suffice 
for  the  determination  of  its  spectral  properties. 
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L  The  Etlect  of  Errors  in  the  Estimate  of  the  Power  Spectral 
Density  Function. 

The  nature  of  the  errors  in  vehicle  parameter  studies  based 
upon  inaccurate  power  spectral  density  forms  and  measurements  is 
of  prime  importance  in  our  studies.  As  one  must  expect,  we  may 
only  obtain  the  spectral  density  function  from  surveys  within  some 
confidence  band. 

In  the  following  paper  two  criteria  are  used  upon  which  to 
make  optimum  judgements  concerning  the  parameters  that  charac¬ 
terize  the  simple  vehicles  considered.  The  criteria  concern  the 
variance  of  the  vertical  acceleration  and  the  stationary  values  of  the 
power  spectral  density  of  the  vertical  acceleration. 

One  can  easily  show  that  the  power  spectral  density  of  the 
output  of  a  linear  system  which  is  subjected  to  a  stationary  random 
input  is  simply,  assuming  one  dimensional  spectra 

26)  P(X)  =  \H(iX)\*  P(X)  =  h(X)  P(X) 

where  H(iX)  is  the  frequency  response  of  the  linear  system. 

We  must  recognize  that  h(\)  is  really  A  ( A  ,/cf )  where  T  •*> 
the  parameter  vet  tor  whose  values  we  wish  to  optimize  relative  lo 
the  criteria  mentioned  above.  We  shall  consider  first  order  otlcets 
on  the  second  of  the  criteria  mentioned  above  from  errors  in  the 
power  pe.  tral  density  function. 

#  0 

taniS'.dei  the  power  spectral  dens.it>  P(\)+P( X),  whr  re  P(X)  is 


20 


the  perturbation  or  error  term.  The  maximum  values  of 
P(X)  [  P(X)  +  rtxj]  are  determined  aa  usual  from  the  equation 

27)  [p(A)+  p*(A)\  +  h(A,p)  +  m 

We  have  specialized  the  parametric  study  to  only  one  param¬ 
eter  p  .  We  write  the  solution  to  ?7)  as  A(p)  +  )P(p)  where  A(p) 
is  the  root  if  there  is  no  error  in  the  power  spectral  4*nsity  estimate. 
The  curve  of  the  maxima  as  a  function  of  P  has  the  form 

it)  HAtpi+ttp),  t)  [p(Mp)*f(t)/+  rtMp)*t(p)i\ 

From  27)  we  may  find  the  first  order  estimate  of  )2tp)  by  substitut¬ 
ing  A(p)+f(P)  into  27),  then  expanding  this  expression  about  X  (p)  . 
Upon  using  the  fact  that 

29)  -jf- (Mp)’ pJ  P(X(p))+h(A(p),  p)  -*fr(X(pJ)-  0  , 


we  find 

30)  X(p)— 


d 

dX 

h(  X  (P)  ,p)  rfX  (p)J 

L 

d£ 

[a  X(p),p)  P  X(P))} 

Proceeding  in  the  same  fashion  with  equation  28),  we  can  determine 
the  first  order  difference  in  the  curve  of  the  maxima  as  being 

3D  h(X(p),p)  P'(X  (pJ). 

By  use  of  the  relations  30)  and  31)  we  can  therefore  estimate  the 
magnitude  of  the  errors  involved  by  approximating  the  power  spec¬ 
tral  density  function  as  exponential. 
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VII,  Ground  Surveys  and  Conclusions 


The  Midwest  Applied  Science  Corporation  has  under  contract 
with  the  Land  Locomotion  Laboratories  undertaken  the  task  of  anal¬ 
yzing  actual  ground  survey  data  to  determine  the  nature  of  real  ground 
spectra,  and  make  decisions  regarding  characterizations  of  the  sur¬ 
face.  To  the  best  of  our  knowledge  this  has  not  been  undertaken  in 
the  U.  S.  A.  before.  The  oceanographers  are  certainly  well  ahead 
in  this  approach. 

Since  there  is  no  readily  available  data  for  our  purpose  of  study¬ 
ing  the  ground  as  a  surface,  we  have  performed  a  survey  on  a  field 
that  has  been,  in  its  history,  plowed,  cultivated  and  grazed.  This 
field  is  located  in  North-central  Indiana,  adjacent  to  Purdue  Univer¬ 
sity,  and  is  part  of  the  Purdue  dairy  farms. 

The  survey  on  this  field  was  taken  along  five  concentric  lines 
that  divide  a  circle  into  equal  sections,  as  shown  in  Figure  4. 

The  lengths  of  the  lines  vary  from  fifteen  hundred  feet  for  line 
A  to  twelve  hundred  feet  for  line  D.  The  elevations  were  taken  at  six 
inch  intervals  along  each  line,  thereby  yielding  a  spectral  resolution 

r 

of  one  cycle  per  foot. 
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1800' 


Figure  4.- 

It  was  decided  to  use  the  lines  in  the  configuration  shown  rather 
than  a  rectangular  grid  on  the  surface  due  to  the  fact  that  we  are  pri¬ 
marily  interested,  lor  the  land  locomotion  problem,  in  the  spectral 
properties  as  a  function  of  direction.  Hence  a  polar  survey  appeared 
more  suitable  than  a  rectangular  survey. 

The  use  of  five  lines  was  dictated  by  the  fact  that  if  we  are  to  try 
to  til  ellipses  to  the  constant  power  contours,  then  this  elliptical  distri¬ 
bution  of  the  spectral  properties  with  direction  would  m  general  be  de* 
termined  by  five  constants.  It  is  expected  that  the  center  of  the  ellipse 
will  be  very  near,  ji  not  -m,  the  lrequetic  y  origin.  For  t  his  n  -ison  we 


have  assumed  only  three  arbitrary  constants  in  the  previous  analysis. 
However,  we  shall  continue  to  use  five  lines  for  the  first  few  surface 
surveys . 

To  obtain  the  power  spectral  densities  a  numerical  analysis  was 
performed  on  the  raw  data,  and  programmed  for  computation  on  a  dig¬ 
ital  computer.  These  are  presented  along  with  their  ninety  per  cent 
confidence  limits  and  their  variances  in  figures  5-8.  The  numerical 
procedure  used  on  the  raw  data  was  first  to  subtract  out  least  squares 
parabolic  base  lines  along  each  contour  in  order  to  delete  the  extreme¬ 
ly  low  frequencies,  for  these  frequencies  do  not  effect  the  trafficability 
of  vehicles.  The  resulting  data  were  then  analyzed  according  to  the 
methods  of  Blackman  and  Tukey  |^9j. 

The  method  of  deleting  the  very  low  frequencies  is  still  open  to 
investigation  due  to  the  possibility  of  introducing  high  frequency  com¬ 
ponents  from  the  points  at  which  the  parabolic  base  lanes  are  connect¬ 
ed.  This  may  account  for  the  large  power  concentrated  around  .5  cycles/ 
ft.  on  line  B,  and  on  the  large  variations  of  the  computed  variances. 
Clearly,  the  base  lines  must  effect  these  variances. 

The  complete  analysis  of  these  power  spectra  is  presently  in 
progress,  therefore  it  is  too  early  to  say  whether  this  data  fits  the 
models  de-scribed  previously.  Preliminary  studies  have  indicated  a 
change  in  the  functional  form  to  fit  the  c  urve  at  the  lower  lrequencios 
more  faithfully.  But  since  this  is  intimately  connected  with  the  .nitial 
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data  redaction  procedures  we  cannot  make  a  decision  on  this  point 
until  further  studies  are  performed. 

It  is  our  feeling  that  the  future  will  see  the  off-road  vehicle 
that  will  be  adaptively  controlled  to  adjust  itself  to  the  terrain  that 
it  must  operate  on.  The  signal  that  the  control  mechanism  must 
operate  on  is  derived  from  those  aspects  of  the  ground  surface  that 
will  be  continuously  measured,  or  sampled,  and  the  optimizing  cri¬ 
terion  that  is  built  into  the  system.  But  what  aspects  of  the  ground 
surface  contour  should  be  measured  remembering  that  it  is  a  ran¬ 
dom  surface?  We  feel  that  this  can  be  answered  by  studies  along 
the  lines  presented  in  this  paper.  We  hasten  to  add  that  this  is  just 
an  initial  study  of  the  general  probelm  from  the  off-road  mobility 
point  of  view.  Therefore  while  particulars  may  change,  the  general 
philosophy  presented  will  remain  fixed. 
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